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Abstract 

We investigate the possibility that the galactic dark matter exists in an scenario where the phan- 
tom field is responsible for the dark energy. We obtain the statically and spherically approximate 
solution for this kind of the galaxy system with a supermassive black hole at its center. The so- 
lution of the metric functions is satisfied with gu = —grr- Constrained by the observation of the 
rotational stars moving in circular orbits with nearly constant tangential speed in a spiral galaxy, 
the background of the phantom field which is spatially inhomogeneous has an exponential potential. 
To avoid the well-known quantum instability of the vacuum at high frequencies, the phantom field 
defined in an effective theory is valid only at low energies. Under this assumption, we further inves- 
tigate the following properties. The absorption cross section of the low-energy 5- wave excitations 
of the phantom field into the central black hole is shown to be the horizontal area of the central 
black hole. Because the infalling phantom particles have a total negative energy, the accretion of 
the phantom energy is related to the decrease of the black hole mass which is estimated to be much 
less than a solar mass in the lifetime of the Universe. Using a simple model with the cold dark 
matter very weakly coupled to the " low-frequency" phantom particles which are generated from the 
background, we show that these two densities can be quasi-stable in the galaxy. 
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I. INTRODUCTION 



The recent experimental data have shown that the current Universe is undergoing a phase of 
accelerated expansion Considering the universe filled with a barotropic perfect fluid which 

corresponds to the dark energy component, its equation of state w < —1/3 is required for cosmic 
acceleration, where w = p/p with p and p being the density and pressure, respectively. Recently 
observations suggest that the equation of state lies in a narrow strip around w = —1 0, Q| j where 
w = —1 corresponds to a cosmological constant A and < — 1 is allowed. A specific form of the 
dark energy corresponding to the phantom field was proposed to realize the possibility of late-time 
acceleration with w < —1 0, @], while the quitom model has w crossing —1 The peculiar 
property of the phantom dark energy is the violation of the dominant energy condition (DEC), so 
that the energy density and curvature may grow to infinity in a finite time, which is referred to a 
Big Rip singularity 0,]^. 

Observations related to the cosmic microwave background (CMB) and the large-scale structure 
(LSS) support that the Universe is very close to spatially flat geometry 0,0. At the present time 
our universe is dominated by dark energy with the fraction ~ 72%. The most accessible component 
of the Universe is baryonic matter which amounts to only 4.6%. The main remaining part that is 
non-baryonic and non-luminous is believed to be the so-called dark matter responsible for ~ 23%. 
The weakly interacting massive particles (WIMPs) are considered one of the main candidates for 
cold dark matter which is dust-like with equation of state w 0. A dark matter halo, demonstrated 
by its gravitational effect on a spiral galaxy's rotation curve, dominates the mass of a galaxy and 
can stretch up to be larger than 50 kpc from the center of a galaxy. 

The presence of the interaction between the dark matter and dark energy field may modify the 
distribution of dark matter in a galaxy. A homogeneous scalar field with negative kinetic energy has 
)een used to model the phantom energy (with w < —1) and then to study the cosmological evolution 



1CH12||. where the scalar field varies only with time but does not change spatially. However, for a 



galaxy due to the gravitational instabilities the phantom field may not be spatially hor nog eneous. 
Some works related to the inhomogeneous dark energy properties have been completed 13l42ll|. In 
this paper, we are interested in the static solution of the Einstein equations that can describe the 
dark matter halo with the existence of the supermassive black hole at its center and the background 
of the spatially inhomogeneous phantom field. We adopt the standard assumption that the dark 
matter halo consists of the WIMPs with ~ 0. We find that an approximate solution of the metric 
exists for describing this galactic halo scenario when we take the metric functions to be satisfied 
with gtt = —Qrr- general, several static and spherically symmetric exact solutions of Einstein 
equations are obtained in terms of the parameter e, defined by T^g = T\{1 — e), for which some 
related works can be found in Refs. 22h25||. We also obtain the approximate solution for the 



spatially inhomogeneous phantom field in a galaxy. Our result indicates that the corresponding 
exponential potential of the phantom field is relevant to the stage of structure formation. 

We further study the stability of the space-time structure for the galaxy. Since the phantom 
dark energy field can cause the quantum instability of the vacuum at high frequencies, we therefore 
treat it as an effective theory valid at low energy, i.e., we add a cutoff in the momentum integral 
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[lo| . The stringent limit for the cutoff is A ^ 3 MeV, which was obtained from the diffuse gamma 
ray background 26|]. We compute the accretion rate of the phantom particles into the black hole. 
Since the infalling phantom particles have a total negative energy, the black hole mass diminishes in 
the process. On the other hand, due to the conservation of the angular momentum of the individual 
WIMPs and their tiny interaction rate, the capture cross section for dark matter particles by the 
supermassive black hole is sufficiently small. Therefore, the dark matter that we consider here is 
quite stable. 

This paper is organized as follows. In Sec. [Ill we consider the pressless massive dark matter 
(WIMP) directly or indirectly interacting with a phantom field which is associated with the ac- 
celeration of the Universe. We solve the Einstein equations with the condition gtt = —Qrr, and 
the approximate solution is then obtained. As a byproduct, we will show that the exact solutions 
can be extended to some typical limits which are related to Schwarzschild, Reissner-Nordstrom 
and Schwarzschild-de Sitter /anti-de Sitter solutions, respectively. We calculate the distribution 
function of the phantom field and its potential in the region of the galactic halo. In Sec. IIIH we 
will examine if the spatially inhomogeneous phantom field can stably survive with the existence 
of supermassive black hole at the center. Considering the spherically symmetric space with a su- 
permassive black hole at the center, we first semiclassically calculate the Klein-Gordon equation of 
a 5-wave phantom particle that is excited from the background field. We obtain the absorption 
probability of the infalling phantom particle into the central black hole and further show that the 
absorptive cross section is approximately the area of the horizon. We show that the accretion rate 
of the phantom particle, which is accompanied by the decreasing rate of the black hole mass could 
be small enough, so that the space structure of a galaxy is stable compared with the life of the 
Universe. In Sec. IIVI we consider a toy model of phantom particles coupled to massive dark matter, 
and show that both the dark matter and low frequency phantom densities can be quasi-stable for 
a sufficiently small coupling constant. Finally we give the summary in Sec. |Vl 



II. THE EXACT SOLUTION IN THE STATIC LIMIT 



We consider the real phantom field minimally coupled to gravity 
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where = 8ttG is the reduced Planck mass, V{^) is the phantom field potential, the Lagrangian 
term Crn accounts for the massive dark matter in the galaxy, and Cj describes possible interac- 
tions between the phantom field and the dark matter. Due to the small coupling between the 
phantom field and the dark matter, Cj can be negligible in the present calculation. To investigate 
static, spherically symmetric solutions, we employ the metric ds"^ = —e^dt^ + e^dr"^ + r'^{dO'^ + 
s\v? 9d<jP'), adding the ansatz \ = —v. This condition is satisfied with the exact solutions, like the 



Schwarzschild, Reissner-Nordstrom and de Sitter /anti-de Sitter solutions. For the static situation 
the Einstein equations read 
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where the energy-momentum tensor corresponds to the massive dark matter in the background of 
the phantom field, 
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and a prime denotes the differentiation with respect to r. In the spherical coordinate, the to- 
tal energy-momentum tensor is denoted by diag(— the phantom field background is 
diag{— Pph, Pr,phT Pe,ph, Pcf>,ph) ^ and the cold dark matter can be approximated as diag(— pdm, 0, 0, 0). 

It is interesting to note that we can find solutions that are satisfied with the condition A = —ly, 
i-e. gtt = —Qrr, in the limit T\ — )• T^^.. Moreover, it will be shown in the following that as for the 
mass density of WIMPs to be pdm = e^^^'^ > 0, the present case can be satisfied with this condition 
for which the existence of the corresponding solution is due to the fact that the sign of the kinetic 
term of the phantom field is opposite compared to the ordinary scalar field (quintessence field); for 
the quintessence dark energy, there exists no such solution. 

To find the solution satisfied with the condition A = —u, we first set u = ln(l — U) and substitute 
it into Eqs. ([2]) and O, (or into ([3]) and (P). We then obtain 



r^U" + 2erU' + 2{e - 1)U = , 
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where we have set T^q = T^, = T*j(l — e) with e being a constant. The solution of this equation is 
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where rs,r^, and a 7^ are the integration constants. 

Before we continue, we discuss the obtained exact solutions in some typical limits. In addition to 
the present case of the galactic dark matter interacting with the phantom field, the following ones 
can be satisfied with the relation = T\. First, for e = 1 with — t- 00, it gives the Schwarzschild 
metric corresponding to T'^ = 0, and is the Schwarzschild radius. Second, for e = 2, the 
solution is the Reissner-Nordstrom metric, where = GQ^ and Q is the charge of the black hole. 
Third, for e = resulting in T*^ = T\ = T^g = T*^ = 3/(87rGre), it gives the Schwarzschild-de 
Sitter/anti-de Sitter solutions which are equivalent to the replacement = — 3/A, with A being 
the positive/negative cosmology constants. 

In this paper, we are interested in the metric that can describe the motions of stars in the galaxy. 
For the several observed cases, the rotational stars with radius rhaio > t ^ in a spiral galaxy, 
where rhaio ~ 50 kpc denotes the radius of a typical halo in a galaxy, are moving in circular orbits 
with nearly constant tangential speed v which roughly ranges from 10~^ to 10^'^. For the region 
with r > rhaio) the dark matter may become very dilute. In the dark matter dominant region, 
where the test particle stably moves in constant rotational curve, the metric function —gu was 
estimated in the form 27 1 
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where rg is a constant and —gu = ^'^ will be denoted as / in the following discussion. Note that 



the form of gu is model independent but grr is not 27|, |28|. Because 2u^ ~ 10 ^ — 10 ^ is a tiny 



magnitude, we can approximate the metric function in the form 
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On the other hand, we focus on our solution with e ^ 1 but 7^ 1, for which the galaxy has a 
supermassive black hole of the Schwarzschild radius r<j ~ 10^^ pc at its center and is predominated 
by the massive dark matter (WIMPs) in the phantom dark energy background. In the following, 
we will further exhibit that if we put a test particle in such a galaxy, it is possible to find a solution 
for which the test particle moves in circular orbit with nearly constant tangential speed consistent 
with the result given in Eq. p6|) (see the result given in Fig. [T]). Compared with Eq. (|16|) . we find 
that in the same dark-matter dominant region if the metric function in Eq. (jl3p is given by 
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where | In (^^^ | < 1, I7I, |a| <C 1, 1 — e = a/2, r^^ = 7/fo , these two equations are approximately 
equal under the following conditions: 7a = 2v'^ ~ 10~^ — 10~^ and 7 = 2v'^ ln(fo/ro). In Fig. [H we 
consider the extreme case of the dark-matter galaxy with the supermassive black hole at its center 
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FIG. 1. The metric function —gu as a function of r. The soUd curve corresponds to —gtt = 
1 — ^ + 7(^)° and the dashed curve is iox —gtt = [■^\ , where we have adopted the fohowing 
parameters: = 10"^ pc, rg = 6 Mpc, a = —0.155, = 10~^. For comparison, using the 
same rg, the upper and lower dot-dashed curves are for —gtt = (^^^ with v"^ = 10-''/4 and 

v"^ = 4 X 10"^, respectively. 

and show the corresponding metric function —gtt as a function of r. We find that the obtained 
metric function can be consistent with the result given in Eq. p5p under the certain condition. 

In the region with the nearly constant tangential speed of the stars, the energy density is given 
by /) ~ ?;^mpj/47rr^, where nipi = G~^^'^ is the Planck mass. Moreover, from Eqs. ([2]), (j3]), ([6]), and 
([7]), we have Pdm = e'^^'^ in the static limit, where is the classical phantom field background. On 
the other hand, from Eqs. ([6]), ([7]) and dH) and using e 1, we find that p = e'^^'^/e ~ />dm = e'^^'b^ 
i.e., pph ~ and V = + p(l — e) ~ Since the metric function is e'^ ~ 1, combining 

these gives 



$,(r)-$oo^^In(^) , (18) 
\/47r \roJ 

, (19) 

during the distances rtaio > r ^ Vg, where the integration constant fo is roughly ^ 6 Mpc, the 
intergalactic distance. Here we choose the positive sign of <!>{,. As for r ^ fo, the space becomes 
flat, the phantom field = ^oo is spatially uniform, and its potential is responsible for the current 
accelerated expansion for which we have S-ffg ~ 87ry(<l>oo)/"ipi, where at the present epoch the 
Hubble parameter Hq ^ 10~^^ GeV. Consistently, we get 

fo ~ ~ 6 Mpc. (20) 
-no 



r 



III. THE ABSORPTION OF THE PHANTOM FIELD BY THE SUPERMASSIVE 
BLACK HOLE 



The WIMPs (dark matter) are quite stable owing to the conservation of the angular momentum 
of the individual particles and the tiny interaction among them, so that the capture cross section 
for dark matter particles by the supermassive black hole is sufficiently small. Therefore, to study 
the stability of the space-time structure for the galaxy, we compute the accretion rate of the excited 
phantom particles into the black hole. If the non-relativistic dark matter can be long lived enough 
compared with the age of the universe, in the following, we will estimate the absorption probability, 
absorption cross section, and the accretion rate of the phantom excitation wave into the black hole. 
At low energies, the dominant effect comes from the excited phantom particle with lowest angular 
momentum. Therefore, we consider its S state, which is the excitation from the background, 
5$ = <&(i, r) — ^*b(r). The phantom potential is assumed to be around the local maximum, 

-pi , (21) 

where -F(<5) — + (-^*$(^6)/2)(J^*^ (since the phantom particle with the negative kinetic energy 
might evolve to the maximum of the local potential). In other words, the local potential is around 
the locally stable point, V$(<I>5) = 0. Here F<j>$ = cPF/d^'^ and V$ = dV/d^. Taking into account 
the back-reaction, the wave for the phantom excitation is satisfied with the Klein-Gordon equation 

- J^oS^ + ^drir^fdrS^) = (ml + X{^^))6^ = mlsS^ , (22) 



where 



2 T^/;^ ^ \/l67r 



^ vnipi 



is the mass squared of the excited phantom field S^, A is the dimensionless coupling constant 
describing the interaction between the excited phantom particle and the dark matter in the form of 
jCi = (A/2)(5$2 . ^q,2^ stands for the average of the quantum fluctuation of the dark matter 

field ^. We will show in the next section that, when A is small enough, the densities pg^ and p>j. 
exhibit a stable oscillatory behavior, so that one can take m^g ~ m^. In this paper, the conclusion 
is also applicable for m| < 0. 

The excited phantom wave can be variable separated in form of S^{t,r) = 3^^[^(r)e~*'*'*]. In 
the present study, we reasonably assume the tiny phantom mass and therefore neglect its spatial 
dependence. We redefine a new coordinate variable r by the relation 

= . (23) 

Then, the differential equation of <p reads 

[d', + {co'-mlsfy]ct> = 0. (24) 



Very close the horizon r ~ Vg, where / ~ and we can have the approximation (a; — mgg/)r' 
Lo'^rg, the ingoing wave is then given by (f){r) ~ (j)Qe~'^^'^^'^ . In this case, because the infalling phantom 
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particles have a total negative energy, the phantom energy accretes resulting in the decrease of the 
291 ] ■ Extending to the region, where r ^ wr^ and the rotational speed, v, of the 



black hole mass 



star is constant, corresponding to / ~ (r/rg) [27[, the solution can be approximated as 
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where d ~ [1 — rrigg/w^J^/^. On the other hand, using the dimensionless variable p = ur, we find 
that the differential equation for (p can be rewritten by 
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For the distances rhaio > r S> r^, this equation is approximated as 
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The corresponding solution is 



(26) 



(27) 



(28) 



The overlap region between the two solutions, Eqs. (j25p and (j28p . is r| ^ <C 1, where Eq, 
(l28l) reduces to 
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Matching Eq. (|29p onto Eq. (j25p in the overlap region, we find 
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r(| + ^2)r(i+„2),-2.(i + 2^2) 

From the above result, we expect \B\ ^ \A\ for ujTs <ti 1. For distances ojr ^ 1, the solution in 
Eq. (j28p asymptotically behaves as 

1 



\/27rd 
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where 6 = dujr — vru /2 — 7r/2. Thus, we know that the absorption probability of a spherical S wave 
is 

2 
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A6 



(32) 



We then further calculate the absorption cross section for the ingoing S wave. For an incident 
plane wave (/>(z) = e*'^^, which can be expanded in the partial- wave amplitudes 



Akr cos 9 



^(2£ + l)P^(cos( 



ikr —i{kr—i-iT) 



(33) 



8 



the absorption cross section of the 5-wave component is 

number of particles absorbed by the area of spherical surface per unit time 



fabs 



number of incident particles crossing unit area per unit time 



2kr 
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(34) 



1 A;2 

where k = duj, r is the area of the spherical surface, and A-ktI is exactly the area of the Schwarzschild 
horizon. The result is consistent with the low-energy cross section for massless minimally coupled 
scalars [30]. It should be stressed that the result is independent of r. Because the phantom field 
has a negative kinetic term, the phantom energy flux onto the black hole is Tq^ = —6^^t^^,r 
which has the opposite sign compared with the ordinary matter fluid. For the ultralight |7n||, the 
potential term can be negligible for r < r^aio, and we can approximately take the Jacobson solution 



(5$ = 6^oo[t + i"s lii(l — Tg/r)] near the horizon 31[, where 5^oo is the excitation of the phantom 
field in the absence of the black hole. The excited phantom particles carry negative energies. 
As they fall into the black hole, the black hole diminishes with the rate dMBu/dt = Anr'^Tor = 
—47rrg(5$oo)^i accompanied by the accretion of the phantom energy 29j, |33, l33| . In 32l. 1331] . the full 
nonlinear absorption of a phantom field by a black hole has been taken into account. Alternatively, 
the above calculation can be obtained using the Eddington-Finkelstein coordinates for which the 
phantom energy flux across the absorption area (horizon) is T^y = — (5$oo)^) where the advance 
time coordinate u = t + r + In ^^^^ . Thus the decrease of the black hole mass with the same rate 

The boundary condition far away the galaxy center is based on the assumption that the phantom 
field is spatially uniform. We adopt the simple potential model V{^) = m?^'^ with m ~ 10^'^'^ 



eV [3j] to describe the acceleration phase at the present epoch. We can match this model with 
Eq. p9|) . when $ is very close to ~ ""^ph If the kinetic energy of the phantom field tends to 
be subdominant compared with its potential energy, from the equation of motion we can have the 
approximation <i>oo — V$/(3if) = mmpi/\/67r. In this model no Big Rip occurs and it is satisfied 
with 1(7 — >• — 1 for t — 7- oo. It is expected that the perturbations 6^oo of the phantom field are of order 
10~^$oo Therefore, if the kinetic energy of the excited phantom states is also expected to be 
much less than or even the same order of magnitude as that of the background, ((5$oo)^ ~ C'($^), 
the decrease rate of the black hole can be estimated to be dMBu/dt ^ — IO^^^Mq yr~^, where we 
have used Mbh ^ W^Mq. Therefore, the decrease of the black hole mass is much less than a solar 
mass in the lifetime of the Universe. 



IV. STABILITY 



In this section, we consider the classical evolution of a simple system for which a nonrelativistic 
scalar dark matter couples to the excitation of the phantom field in a galaxy. Note that to avoid 
the quantum instability of the vacuum at high frequencies, the phantom dark energy field defined 



in an effective theory is valid at low energy [10|, l26|] . The stringent limit for the momentum cutoff. 



which was estimated from the diffuse gamma ray background, is less than 3 MeV 2a]. As for 
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the phantom field defined in an effective theory at low energies, we will exhibit that this system 
can be quasi-stable for a weak coupling. It should be noted that even for low frequencies, the 
negative energy of phantom particles may cause the system to become unstable since the positive 
energy dark matter particles could increase the energy to any level as long as the phantom particles 
decrease the same magnitude of the energy. A similar &TOtem with the phantom excitation coupled 
to the massless graviton had been considered in Ref. 10|. 
The relevant Lagrangian C, is 



where ^ is the nonrelativisitic dark matter field and 5$ is the excitation of the phantom field as in 
the previous study. The spatial variations of them are relatively small in this consideration because 
^ is nonrelativisitc and the accretion rate of the phantom excitations is typically small. Thus the 
energy densities are 



(36) 



and the interaction term /Oiji^^ = — (A/2)^'2(5$2. Note that here summarize only for low-energy 
modes. The equations of motion of these two kinds of particles are 



1 - 

2 



ra. 



5$, 



(37) 



where we have used the dimensionless variables 



A = \{Mlm^y 



(38) 



where pq, ~ v'^mpy/Anr'^, and in the dark- matter dominant region of the galaxy, we have adopted the 
approximation for the metric function / — 1. The perturbations, 5^ and ^, are expected to be of 
order M ~ 10~^mpi [lo| . Note that in this section a prime denotes the differentiation with respect 
to r. We plot the evolutions of densities in Fig. [2l where we have adopted the initial conditions, 
^' = S~^' = and = 5$ = 1 at r = 0. Although we have used = 10"^^ eV and = 10^^ 
eV as inputs, we note that evolution curves do not obviously change if the condition <^ m-q/ 
is satisfied. We obtain that for < A < 0.37, the densities oscillate with a stable behavior, where 
Ps^ ~ 0. However, for A < 0, /Oij- and —ps^ grow quickly. In short, we could conclude that the 
stability for the galactic cold dark matter existing in the phantom field background is possible. 
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FIG. 2. Evolution of the energy densities of a coupled pair of the dark matter denoted by the solid 
curve and phantom excitation denoted by the dashed curve. The initial condition is ^' = 6^' = 
and ^ = (5$ = 1 starting at r = 0. The energy densities are displaced in units of m\M'^ /2 and 
satisfied with the constraint ps^ + p^^ + p^s^ = constant. 

V. SUMMARY 

We have studied the possibility that the galactic dark matter exists in an scenario where the 
phantom field, responsible for the dark energy, may not be spatially homogeneous in a galaxy. We 
have obtained the statically and spherically approximate solution for this kind of galaxy system 
with a supermassive black hole at its center. The static exact solution of the metric functions is 
satisfied with gu = —Qrr, of which the relation is also consistent with the black hole solutions in 
the vacuum, electromagnetic and cosmological constant sources, corresponding to Schwarzschild, 
Reissner-Nordstrom and Schwarzschild-de Sitter /anti-de Sitter metrics, respectively. 

Constrained by the observation that the rotational stars in a spiral galaxy are moving in circular 
orbits with nearly constant tangential speed, we have obtained that the background of the phantom 
field which is spatially inhomogeneous has (i) the small density pp^ ~ with negligible fluctuation if 
the coupling between the excitation phantom field and dark matter is small, and (ii) an exponential 
potential. 

To avoid the well-known quantum instability of the vacuum at high frequencies, the phantom 
field that we consider here is an effective theory valid at low energies. Under this condition, we 
have computed the absorption cross section of the S'-wave excitations, arising from the phantom 
background, into the central black hole and shown that it is equal to 47rr^, the horizontal area of 
the central black hole. Because the infalling phantom particles have a total negative energy, we 
estimate that the black hole mass thus diminishes at the quite small rate dM-Qn/dt ^ —10~'^^Mq 
yr~^, so that the decrease of the black hole mass is much less than a solar mass in the lifetime of 
the Universe. Furthermore, using a simple model with the cold dark matter very weakly coupled to 
the " low- frequency" phantom particles which are generated from the background, we have shown 
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that the dark matter and phantom densities can be quasi-stable. 
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